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Different mechanisms operate in various regions of the MSSM parameter space to bring the relic density of 
the lightest neutralino, y?, assumed here to be the LSP and thus the Dark Matter (DM) particle, into the range 
allowed by astrophysics and cosmology. These mechanisms include coannihilation with some nearly-degenerate 
next-to-lightest supersymmetric particle (NLSP) such as the lighter stau fi, stop t\ or chargino y^, resonant 
annihilation via direct-channel heavy Higgs bosons H/A, the light Higgs boson h or the Z boson, and enhanced 
annihilation via a larger Higgsino component of the LSP in the focus-point region. These mechanisms typically 
select lower-dimensional subspaces in MSSM scenarios such as the CMSSM, NUHM1, NUHM2 and pMSSMIO. We 
analyze how future LHC and direct DM searches can complement each other in the exploration of the different DM 
mechanisms within these scenarios. We find that the f\ coannihilation regions of the CMSSM, NUHM1, NUHM2 
can largely be explored at the LHC via searches for I/t events and long-lived charged particles, whereas their 
H/A funnel, focus-point and y^ coannihilation regions can largely be explored by the LZ and Darwin DM direct 
detection experiments. We find that the dominant DM mechanism in our pMSSMIO analysis is y* coannihilation: 
parts of its parameter space can be explored by the LHC, and a larger portion by future direct DM searches. 
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l. Introduction 

The density of cold dark matter (CDM) in the 
Universe is now very tightly constrained, in par¬ 
ticular by measurements of the cosmic microwave 
background radiation, which yield Ucdm^ 2 = 
0.1186 ± 0.0020 1 and are consistent with other, 
less precise, determinations. This determination 
of the CDM density at the percent level imposes 
a corresponding constraint on the parameters of 
any model that provides the dominant fraction of 
the CDM density. This is, in particular, true for 
supersymmetric (SUSY) models with conserved 
l?-parity in which the CDM is provided by the 
stable lightest SUSY particle (LSP) [2]. In a 
series of recent papers incorporating the data 
from LHC Run 1 and elsewhere, we have im¬ 
plemented the dark matter (DM) density con¬ 
straint in global analyses of the parameter spaces 
of different variants of the minimal SUSY exten¬ 
sion of the Standard Model (MSSM), assuming 
that the LSP is the lightest neutralino Xi- The 
models studied included the constrained MSSM 
(CMSSM) with universal soft SUSY-breaking pa¬ 
rameters (mo, rap 2 and A q , in standard notation) 
at the GUT scale [3], the NUHM1(2) in which 
universality is relaxed for both together (each sep¬ 
arately) of the soft SUSY-breaking contributions 
to the masses-squared of the Higgs multiplets 
m 2 Hi 2 j3j|4 , and a version of the pMSSMIO |(5j, 
in which TO of the effective Lagrangian parame¬ 
ters (3 gaugino masses Mp 2 , 3 , 2 squark masses 

m. q 12 ^ rriq 3 , a common slepton mass mjr, a com¬ 
mon trilinear coupling Aq, the Higgs mixing pa¬ 
rameter /u, the pseudoscalar Higgs mass Ma, and 
the ratio of Higgs vevs tan /3) are treated as inde¬ 
pendent inputs specified at the electroweak scale. 

Reproducing correctly the cosmological CDM 
density requires, in general, some special choice of 
the SUSY model parameters, which may be some 
particular combination of sparticle masses and/or 
couplings. Examples of the former include hy¬ 
persurfaces in the SUSY parameter space where 
the LSP is almost degenerate in mass with some 
next-to-lightest SUSY particle (NLSP), such as 
the lighter stau f\ 6,7|, stop t\ [§ or chargino 
xf 9 j , or where m^o is almost half the mass of a 
boson such as a heavy Higgs H/A |10], a light 


Higgs h or Z 11 , in which case rapid direct- 


channel annihilation may bring the CDM den¬ 
sity into the allowed range. Examples of special 
coupling combinations include the focus-point re¬ 
gion 12 , where the LSP acquires a significant 
Higgsino component. 

We have commented in our previous work on 
the relevances of these DM mechanisms for our 
global analyses. Here we discuss systematically 
which DM mechanisms are dominant in which 
subspaces of the CMSSM 


13 , NUHM1 14 15 


NUHM2 15, 16 and pMSSMIO 17 parameter 


spaces, what are the corresponding experimental 
signatures, and how one might discover SUSY in 
each of these different DM regions. 

Our analysis of the possible detectability of su¬ 
persymmetry in the CMSSM, NUHM1, NUHM2 
and pMSSMIO, depending on the dominant DM 
mechanisms, is summarized in Table |Tj 


2. Measures of Mass Degeneracy 


We first introduce measures on the MSSM 
parameters that quantify the relevant mass 
degeneracies and define each of the above- 
mentioned subspaces in the CMSSM, NUHM1 


and NUHM2 [§[18 
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In each case we also indicate the colour cod¬ 
ing we use in the subsequent figures. The mea¬ 
sures 0 that we use are empirical, but we have 
verified extensively that CMSSM, NUHM1 and 
NUHM2 points that satisfy the DM density con- 
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straint do fulfill at least one of these conditions, 
and that they indeed correspond to the dominant 
DM mechanisms (in the sense of giving the largest 
fractions of final states, generally > 50%) (4j 18 


We have found that there are some ‘hybrid’ re¬ 
gions where the dominant mechanism requires 
two of these conditions simultaneously. In partic¬ 
ular, there are regions where the dominant DM 


mechanism is 


~+~« 


—> bb or tt, processes involv¬ 


ing both stau coannihilation and annihilation via 
the A/H funnel, which we colour purple. There 
are also regions where the chargino coannihilation 
condition is satisfied simultaneously with the stau 
coannihilation or A/H funnel condition. How¬ 


ever, a dedicated study using MicrOMEGAs 19 


shows that chargino coannihilation is the dom¬ 
inant DM mechanism in these regions, and that 
hybrid processes dependent on the m -± and some 

Xi 


other degeneracy conditions being valid simulta¬ 
neously are unimportant, so we colour these re¬ 
gions the same green as the other chargino coan¬ 
nihilation regions. 

The above DM mechanism conditions need to 
be modified for our analysis of the pMSSMIO. 
First, as we see later, funnels due to annihila¬ 
tions via direct-channel h and Z resonances can 
be important [ll|, so for the pMSSMIO we add 
to (|T]) the supplementary criteria: 


h funnel (magenta) : 


M h 
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< 0.4, 


Z funnel (orange) : 


Mz 


m Xi 



< 0.4. 
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Secondly, we find that chargino coannihilation 
dominates in the pMSSMIO also when the sec¬ 
ond condition in 0 is relaxed: we use later the 
condition \m-±/m^o — 1| < 0.25 in our subse- 
quent analysis, which reproduces better the do¬ 
mains of dominance by y^ coannihilation 0 Fi¬ 
nally, we recall that the focusing property of the 
RGEs is not relevant in the pMSSMIO. However, 
the LSP annihilation rate may still be enhanced 
when the fifth measure in 0 is satisfied, due to 

lr This approach yields results similar to those of |l8] , where 
empirical constraints combining the masses and neutralino 
mixing matrix elements are used. 


a larger Higgsino component in the LSP, though 
we find that the dominant DM mechanism in the 
pMSSMIO generally does not involve this prop¬ 
erty. dWe use the same cyan colour to identify 
regions where this condition is satisfied, though 
it is not due to focus-point behaviour. 

Our discussion here of DM mechanisms is 
based on our previously-published global likeli¬ 
hood analyses of the CMSSM and NUHM1 j3], 
the NUHM2 [I] and the pMSSMIO (5]. The 
reader wishing to know details of our treat¬ 
ments of the various experimental, phenomeno¬ 
logical, theoretical and cosmological constraints, 
as well as our strategies for scanning the parame¬ 
ter spaces of these models is referred to [3}|5l[20 


3. Dominant Dark Matter Mechanisms 

In this section we discuss the various mecha¬ 
nisms that play dominant roles in bringing the 
relic density into the experimentally measured in¬ 
terval in our four models. We display in Fig. [T] 
(mo, m 1 / 2 ) planes for the CMSSM (upper left) |3|7 
the NUHM1 (upper right) [3 and the NUHM2 
(lower left) [4], while for the pMSSMIO we show 
the (mq,m^o) plane (lower right), where to, de¬ 
notes the mass of the squarks of the first two 
generations, which we assume to be common [^] 
The Ay 2 3 = 2.30 and 5.99 contours that we found 
in global fits to these models, corresponding ap¬ 
proximately to the 68 and 95% CL contours, are 
shown as solid red and blue lines, respectively. 
Here and elsewhere, the green stars indicate the 
best-fit points, whose exact locations in some pa¬ 
rameter planes are poorly determined and do not 
carry much useful information, in general, as the 
y 2 minima are quite shallow. Also shown, as 
solid purple lines, is the current 95% CL CMSSM 
exclusion from l/J^ searches at the LHC [^J The 
dashed purple contours in the CMSSM, NUHM1 
and NUHM2 cases show the prospective 5 -a dis- 

2 In some cases, these and subsequent figures may include 
small updates from the versions shown previously [3jj5] , as 
they incorporate the latest implementations of the exper¬ 
imental constraints. 

3 As discussed in |s|4], this exclusion curve can be applied 
to the NUHM1 and NUHM2, also in the range of mo < 0 
shown in these plots, where we interpret negative mo = 
Sign(mg)^/mj). 
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covery reaches for ]£t searches at the LHC with 
3000/fb at 14 TeV, corresponding approximately 
to the 95% CL exclusion sensitivity with 300/fb at 
14 TeV. In the pMSSMIO case the dashed purple 
contour shows the 95% CL exclusion sensitivity 
of the LHC with 3000/fb assuming rrig rrig, 
and the dash-dotted lines bound the correspond¬ 
ing sensitivity region assuming rrig = 4.5 TeV. 
All the planes are colour-coded as listed in ([l]|2|, 
with regions where none of these processes are 
dominant left uncoloured. 

We see in the upper left panel of Fig. [l] that 
three DM mechanisms dominate in the CMSSM: 
t i coannihilation at low nig < 2000 GeV, the 
H/A funnel at larger mo and TO 1 / 2 , and the fo¬ 
cus point at larger mo and smaller m\n■ There 
is also a hybrid fi/A/H region extending up to 
(mo,mi/ 2 ) ~ (2000,2500) GeV. In the case of 
the NUHM1 shown in the upper right panel of 
Fig- 0 there is an analogous fi coannihilation 
region. However, it has a much larger hybrid 
ti/A/H region, which has an extension to low 
m 1/2 ~ 1000 GeV for m 0 < 3000 GeV. On 
the other hand, xy coannihilation dominates in 
a large region with my 2 > 2500 GeV. Here 
/i -C Mi, so that the x°,x 2 and Xi are nearly 
degenerate in mass and the x? has mainly a 
Higgsino composition. A similar Xi coannihila¬ 
tion region is visible in the NUHM2 in the lower 
left panel of Fig. [T] where we also see a more 
extensive fi coannihilation region extending to 
?77o ~ 2000 GeV |/J with a relatively small hybrid 
fi/A/H region. This is the only case where we 
see a region of dominance by t\ coannihilation, in 
islands around (mo,rai/ 2 ) ~ (2000,500) GeV. 

In contrast, as shown in the lower right 
panel of Fig. [T] we found in our version of the 
pMSSMIO [5] that the dominant DM mechanism 
is usually Xi coannihilation, this time with a 


4 Analyses with micrOMEGAs confirm that x 1 coannihila¬ 
tion processes contribute over 50% of the final states in 
most of the green shaded regions for the NUHM1 and 
NUHM2, and generally over 75% for m 1 /2 > 3000 GeV. 
In the NUHM1 and NUHM2 for m ^/2 < 1500 GeV, the fi 
coannihilation criterion is also satisfied (and that region 
was colored hybrid in but fi coannihilation does not 
contribute significantly to the relic density calculation, so 
here we colour it green. 


Bino-like LSP and M\ ~ M 2 . [/] We also note 
in Fig. |T] bands at low 0 where rapid annihi¬ 
lations via the h and Z funnels are dominant. 
Not shown in Fig. [T]are scatterings of points with 

0 > 300 GeV where fi coannihilation can also 

.1 

be important, and of points with ra^o <150 GeV 
where the fifth condition in (|T|) comes into play [^] 
We see later that this condition and the f\ coan¬ 
nihilation mechanism dominate in specific regions 
of other projections of the pMSSMIO parame¬ 
ter space, as does ‘bulk’ x^Xi annihilation where 
none of the conditions ([l] [2]) are satisfied. 

One may also consider the possibility that the 
LSP provides only a fraction of the CDM. A com¬ 
plete discussion of this possibility is beyond the 
scope of this paper, but we note that in some re¬ 
gions, e.g., those dominated by fi or t\ coannihi¬ 
lation, lowering the CDM density requires reduc¬ 
ing the NLSP-LSP mass difference. This would 
also have the effect of reducing correspondingly 
the maximal LSP mass, which would favour spar- 
ticle detection at the LHC. On the other hand, 
direct detection would be more difficult if only a 
small fraction of the galactic halos were composed 


of LSPs. For a recent discussion, see 22 


4. The LHC Sensitivity 

In this section we discuss the prospective 
reaches of future LHC searches and their impacts 
in the contexts of the various preferred DM mech¬ 
anisms. We will see that in many cases the pre¬ 
ferred DM mechanism can be directly probed via 
the appropriate LHC searches, as summarized in 
Table [Q 


4.1. Ifix Searches 

Looking now at the physics reach of the LHC 
wither searches 23 , we see that in the CMSSM 


5 We note, however, that this is the result of an interplay of 
various constraints, in particular the pMSSMIO assump¬ 
tion of universal slept on masses, as discussed in detail in 
Ref. [ 5 ]. A different selection of pMSSM parameters could 
favour regions with different dominant DM mechanisms, 
in general. However, this is left for future analysis, and 
here we simply take over the results of Ref. 15]. 

6 Another micrOMEGAs analysis shows that coannihila¬ 
tion processes contribute over 50% of the final states in 
most of the green shaded region, and generally over 75% 
for m ~o < 250 GeV. 

Xl 
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★ - - CMSSM: best fit, lo-, 2a ★ - - NUHM1: best fit, la, 2a 
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Figure 1. 

The (mo,m 1/2 ) 

planes in 

the CMSSM (upper left), the NUHM1 (upper 

right) and the NUHM2 


(lower left), and the (niq,m^ o) plane in the pMSSMIO. Regions in which different mechanisms bring 
the CDM density into the allowed range are shaded as described in the legend and discussed in the text. 
The red and blue contours are the A\ 2 = 2.30 and 5.99 contours found in global fits to these models, 
corresponding approximately to the 68 and 95% CL contours, with the green stars indicating the best-fit 
points, and the solid purple contours show the current LHC 95% exclusions from Ifir searches. In the 
CMSSM, NUHM1 and NUHM2 cases, the dashed purple contours show the prospective 5-a discovery 
reaches for I/t searches at the LHC with 3000/fb at 14 TeV, corresponding approximately to the 95% 
CL exclusion sensitivity with 300/fb at 14 TeV. In the pMSSMIO case, the dashed purple contour shows 
the 95% CL exclusion sensitivity of the LHC with 3000/fb assuming my my, and the dash-dotted lines 
bound the corresponding sensitivity region assuming my = 4.5 TeV. 
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Table 1 

Summary of the detectability of supersymmetry in the CMSSM, NUHM1, NUHM2 and pMSSMIO models 
at the LHC in searches for If t events, long-lived charged particles (LL) and heavy A/H Higgs bosons, and 
in direct DM search experiments, depending on the dominant mechanism for bringing the DM density 
into the cosmological range. The symbols /, (/) and x indicate good prospects, interesting possibilities 
and poorer prospects, respectively. The symbol - indicates that a DM mechanism is not important for the 
corresponding model. The LHC information is drawn largely from Figs. 00 and [/J and the direct DM 
search information from Fig. [5j 


the preferred f\ stau coannihilation region lies 
just outside the current LHC 95% CL exclusion 
region (solid purple line). On the other hand, the 
A/if-funnel region allowed in the CMSSM at the 
95% CL lies well outside this region. However, 
this is no longer the case in the NUHM1 and 
particularly the NUHM2, where portions of the 
A/if-funnel region lie much closer to the LHC 
95% CL exclusion. This is possible mainly be¬ 
cause the Higgs mass constraint is less restrictive 
in these models. In the CMSSM, the future LHC 
sensitivity estimated in Fig. 1 of [24] (dashed pur¬ 
ple line) covers the region where f\ coannihilation 
is dominant, and a part of the hybrid region. It 
also covers slices of the H/A funnel region and of 
the focus-point region. Similar features are seen 
in the NUHM1, except that the focus-point re¬ 
gion is not visible in this case, but the LHC L/t 
search has no sensitivity in the x^ coannihilation 
region. In the case of the NUHM2, the 1 /t search 
is sensitive to only part of the f\ coannihilation, 
and none of the Xi coannihilation region, but it 
does cover part of the H/A funnel region and all 


the t\ coannihilation region. In the pMSSMIO 
case, parts of the Xi coannihilation region and 
the low-m^o band are accessible to future LHC 

_xi 

searches [j though the future LHC Ifir searches 
are less sensitive if nig m,q (dashed line) than if 
mg = 4.5 TeV (dotted line). Table [l] summarizes 
the observability of Lfr events that we estimate 
in the different models considered, depending on 
the dominant DM mechanism in each case. 

4.2. The Possibility of a Long-Lived 
Charged Sparticle 

In some circumstances, a charged sparticle such 
as the Ti or the %i may be the NLSP, and be only 
slightly more massive than the LSP, so that it can 

7 In principle, these bands are accessible to precise searches 
for invisible Z and h decays. However, we found in a sur¬ 
vey of the pMSSMIO o) plane in Fig. |T| that the 

Z —»■ XiXi branching ratio exceeds the current experi¬ 
mental uncertainty of 1.5 MeV for only a handful of the 
lowest-x 2 points. Likewise, the h —> x?X? branching ra¬ 
tio exceeds 0.1 also for just a handful of points. Invisi¬ 
ble Z and h decays both present opportunities for future 
searches. 
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in principle be long-lived. As we now show, the 
presence of a long-lived fi is indeed a distinctive 
prospective signature in the CMSSM, NUHM1 
and NUHM2. 

Fig. H displays on log-linear scales the one¬ 
dimensional Ay 2 profile likelihoods as functions 
of the ti — X? mass differences in the fi coanni¬ 
hilation regions in the CMSSM (upper left), the 
NUHM1 (upper right) and the NUHM2 (lower 
left). Also shown (in the lower right panel) is the 
corresponding distribution in the pMSSMIO, al¬ 
though in this case the ti coannihilation region 
is disfavoured: in this case we see that the likeli¬ 
hood function increases sharply for o < 

10 GeV, rising to Ay 2 ~ 8 for very small mass 
differences. 

In the CMSSM, NUHM1 and NUHM2 panels 
the one-dimensional y 2 profile likelihood function 
is quite a flat function of m y x — to^o, and there 


are points with Ay 2 < 1 that have very small 
values of this mass difference ~ MeV. Hence, 
it is possible that the f\ might live long enough 
(■Tfi > 400 ns) to appear at the LHC as a long- 
lived (LL) charged particle, or even long enough 
(Tfj > 1000 s) to affect Big Bang nucleosyn¬ 


thesis 25-27 . However, we emphasize that the 


mass differences required to realize these possibil- 
ities (< 1.2 GeV, < 0.1 GeV) require quite spe¬ 
cial parameter sets. One should presumably re¬ 
quire ry < 1000 s in order to avoid destroying the 
success of Big Bang nucleosynthesis, a constraint 
that we impose in the following Figures. 


If 


— nr.ro < 1.2 GeV, corresponding to 


7>i > 400 ns, the t-\ would live long enough to be 
detectable at the LHC as a LL charged particle 
29 . Fig. [3] displays the regions of the ( 7770 , 7711 / 2 ) 
plane in the CMSSM (upper left panel), the 
NUHM1 (upper right panel) and the NUHM2 
(lower left panel) where the lowest-y 2 point has 
10 3 s > t y > 10 _7 s: the lifetimes of these points 
are colour-coded, as indicated in the legends. The 
contours for Ay 2 = 2.30(5.99) relative to the ab¬ 
solute minimum of our data set are shown as solid 
red and blue lines, respectively. 

On the other hand, the lower right panel of 
Fig. i displays the one-dimensional Ay 2 func¬ 
tion in the pMSSMIO for the lifetime of the f\ 


in the range 10 3 s > ry > 10 _7 s. We see that 
Ay 2 > 6 throughout the displayed range, in¬ 
dicating that a long-lived fi is not expected in 
the region of the pMSSMIO parameter space that 
is favoured by present data. This is because in 
our analysis the (g — 2)^ measurement and the 
DM constraint favour light sfermions and hig- 
gsinos, whereas long-lived charginos typically re¬ 
quire sfermion and higgsino masses larger than 10 


TeV 


28 


0 


The sensitivity of the LHC to a long-lived fi 
has been compared in 29 to the sensitivity to 
]£t events, and found to be comparable within 
the uncertainties. We therefore assume that the 
projected sensitivity of the LHC to ]pr events is 
a good approximation to its sensitivity to pa¬ 
rameter sets in the fi coannihilation region with 
0.1 GeV < my, — m^o < 1.2 GeV. However, 
while the reach in the parameter space is simi¬ 
lar, the long-lived stau would constitute a spec¬ 
tacular additional signature, and would give com¬ 
plementary information to the direct searches for 
colored sparticles. The purple contours in the 
CMSSM, NUHM1 and NUHM2 panels of Fig. [3] 
again show the present and prospective reaches 
of the LHC for such events. We infer that long- 
lived charged particles could be detectable at the 
LHC throughout the lower parts of the allowed 
regions in the CMSSM, NUHM1 and NUHM2. 
Table [T| also summarizes the observability of long- 
lived charged sparticles in the different models 
considered. 


4.3. Squark and Gluino Searches 

Fig. [4] displays the {m^rrig) planes for the 
CMSSM (upper left), the NUHM1 (upper right), 
the NUHM2 (lower left) and the pMSSMIO (lower 
right). In each panel we show the Ay 2 = 2.30 
and 5.99 contours as red and blue solid lines, re¬ 
spectively. The current 95% CL exclusions from 
ATLAS IpT searches are shown as solid purple 
lines, and the estimated reaches of^r searches for 


95% exclusion with 300/fb of data at 14 TeV 30 


(very similar to the reaches for 5-cr discovery with 


8 The Ax 2 for a X\ in our pMSSMIO sample to have a life¬ 
time in the range 10 3 s > t_± > 10 "'s is much greater, so 
Xi 

this is an even weaker candidate to be a long-lived charged 
NLSP. 










Figure 2. The one-dimensional Ay 2 profile likelihood functions in the CMSSM (upper left), the NUHM1 
(upper right), the NUHM2 (lower left) and the pMSSMIO (lower right) for — m^o < 100 GeV. In the 
CMSSM, NUHM1 and NUHM2, low values of x 2 are found for points with — m^o ~ MeV, whereas 
in the pMSSMIO Ay 2 rises to ~ 8 at small o. 


3000/fb) are shown as dashed purple lines. The 
CMSSM panel shows again that the fi coannihila¬ 
tion region is within the LHC reach in this model. 
However, in the NUHM1 and the NUHM2 only 
portions of the fi coannihilation regions are ac¬ 
cessible at the LHC, along with small pieces of the 
H/A funnel regions. In the case of the NUHM2 
the small t\ coannihilation regions are also well 
within the LHC reach. 

The pMSSMIO panel shows a completely dif¬ 


ferent picture: Xi coannihilation dominates 
throughout the (m,q,mg) plane, as discussed at 
the end of Sect. [31 and the likelihood function is 
very flat beyond the current LHC Ifr exclusion. 
The LHC at 14 TeV will explore a large part of 
the ( mq,mg ) plane, but a (more) complete explo¬ 
ration would be a task for a higher-energy col¬ 
lider [24]. 
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Figure 3. The (too,TOi/ 2 ) planes in the CMSSM (upper left), the NUHM1 (upper right) and the NUHM2 
(lower left), showing the regions where the lowest~x 2 points have Ax 2 < 5.99 and 10 3 s > Tf x > 10~ 7 s: 
the lifetimes of these points are colour-coded, as indicated in the legends. The red and blue contours 
are for Ax 2 < 2.30(5.99) relative to the absolute minimum. Also shown in these panels as solid purple 
contours are the current LHC 95% exclusions from If t searches in the fi coannihilation regions, and as 
dashed purple contours the prospective 5-cr discovery reaches for If t searches at the LHC with 3000/fb 
at 14 TeV, corresponding approximately to the 95% CL exclusion sensitivity with 300/fb at 14 TeV. As 
discussed in the text, the sensitivities of LHC searches for metastable fi ’s in the fi coannihilation region 
are expected to be similar \2fAj . The lower right panel shows the one-dimensional Ax 2 function in the 
pMSSMIO for the lifetime of the fi in the range 10 3 s > tnlsp > 10 _7 s. 


4.4. Charginos and Neutralinos 

The differences between the dominant DM 
mechanisms in the pMSSMIO and the other mod¬ 
els studied are highlighted in Fig. [5j which dis¬ 
plays the (m-± , m^o) planes in the CMSSM (up¬ 
per left), the NUHM1 (upper right), the NUHM2 
(lower left) and the pMSSMIO (lower right). The 


diagonal dashed lines indicate where rn^± = m^o. 
As usual, the red and blue solid lines are the 
Ay 2 = 2.30 and 5.99 contours. In the pMSSMIO 
case, the region preferred at the 68% CL is a nar¬ 
row strip where m-± — m^o is small, whereas in 
r Xi *1 

the other models much of the 68% CL region is 
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Figure 4. The ( mq,rrig ) planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 
(lower left) and the pMSSMIO (lower right). The red and blue solid lines are the Ay 2 = 2.30 and 5.99 
contours, and the solid (dashed) purple lines are the current and (projected) 95% exclusion contours for 
IpT searches at the LHC (with 300/fb of data at If TeV). The solid lines are almost identical with the 
contours for 5-a discovery with 3000/fb. 


in a narrow strip where m-± ~ 2m^o. This re¬ 
flects the fact that in the CMSSM, NUHM1 and 
NUHM2 universal boundary conditions are im¬ 
posed on the gaugino masses at the GUT scale. 

We see that ft coannihilation dominates over 
most of the 95% CL region in this projection of 
the pMSSMIO parameter space, though not in 
the 68% CL region, which has small to-± — m^o 


and is where xf coannihilation dominates [^] On 
the other hand, in the CMSSM, NUHM1 and 
NUHM2, the H/A funnel dominates most of the 
95% CL regions in the (m-±,m^o) planes, and 
also the 68% CL region in the CMSSM, whereas 

®We also note in Fig. [ 5 ] the appearance of an uncoloured 
region with m ~o < 150 GeV, which is where ‘bulk’ anni- 
hilation dominates. 
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Ti coannihilation dominates the 68% CL region 
in the NUHM1 and the NUHM2, as shown in 
Fig. □ At the 95% CL there are also small xf 
coannihilation regions in the CMSSM, NUHM1 
and NUHM2 where m-± — m^o is small, and the 
CMSSM and NUHM1 also have small focus-point 
regions. 

Fig. [5] also displays the present and prospective 
future sensitivities of various LHC searches in the 
(m~±,TO^o) planes. The solid (dashed) orange 
lines are the current and projected 3000/fb 95% 
CL exclusion sensitivities for X1X2 ~* W/Z + 1 / T 
searches)^] , the green dashed lines show the 
projected 3000/fb 95% CL exclusion sensitivities 
for the Xi X2 W/h +$t search, both taken 
from 12l], and the magenta dashed line is the 
projected 3000/fb 95% CL exclusion sensitivity 
f° r X 1 X 21 X 1 X 1 -> T , t -t 2, 3r's +I£t searches, 
taken from 42 . We see that these searches have 


very limited sensitivities to the CMSSM, NUHM1 
and NUHM2, but could explore significant parts 
of the fi coannihilation region in the pMSSMIO. 
However, they would largely miss the xf coanni¬ 
hilation region, i.e., the dominant pMSSMIO DM 
mechanism cannot be explored by direct searches 
at the LHC. 


4.5. The Lighter Stop Squark 

Next we study the differences in the impacts of 
the dominant DM mechanisms on the pMSSMIO 
and the other models in the (m^, m^o) planes 
shown in Fig. [6j In each of these planes, we indi¬ 
cate by dashed lines where = m^o and where 
m= mt + jn^o. In the CMSSM, as shown in 
the upper left panel, we see that the f\ coanni¬ 
hilation region (which contains all of the param¬ 
eter space that is allowed at the 68% CL) is well 
separated from the H/A funnel region, and that 
only a small part of the displayed portion of the 
(m^, m^o ) plane is in the hybrid region. In this 
model we find that — m^o > 300 GeV at the 
95% CL, and we do not find a ti coannihilation 
region, but we do see a focus-point region and a 
small Xi coannihilation region. The situation in 


10 These sensitivities assume that the jq and y) decay ex¬ 
clusively into the x$ in association with W and Z, respec¬ 
tively, not taking into account the decay X 2 X'% 


40 41 


the NUHM1 (upper right panel of Fig. [6]) exhibits 
significant differences. The f\ coannihilation re¬ 
gion (which again dominates the 68% CL region) 
and the H/A funnel region still dominate the dis¬ 
played portion of the (m^ , 0 ) plane, but there 
is a larger hybrid region, the focus-point region 
has disappeared and the xf coannihilation region 
has remained small, but has moved to larger 0 . 
We also note the appearance of a small t\ coan¬ 
nihilation ‘island’ at the 95% CL in this model. 
In the case of the NUHM2 (lower left panel), the 
68% CL region is dominated by f\ coannihila¬ 
tion, whereas its extension to the 95% CL is dom¬ 
inated by the H/A funnel, with small areas of xf 
coannihilation. In this case there is a much more 
prominent t\ coannihilation strip at the 95% CL. 
Finally, in the pMSSMIO (lower right panel), we 
see two patches with small mq — ° connected 

by a narrow ‘isthmus’, and a ‘continental’ region 
at large , where the 68% CL is located. As 
indicated by the green and pink shadings, the 
dominant DM mechanisms in the ‘islands’ are xf 
and fi coannihilation, rather than t\ coannihila¬ 
tion. We also note the reappearance of the h and 
Z funnel bands at low m^o. In this model, the 

Xi _ ’ 

lower-mass ‘island’ and part of the higher-mass 
‘island’ can be explored by future LHC searches 
for h —t bxf jHj, since m-± ~ m^o. This search 
channel is less powerful in the CMSSM, NUHM1 
and NUHM2, where m^± > 2m^o in general, par¬ 
ticularly in the 74 coannihilation regions that are 
favoured at the 68% CL, as seen in Fig. [ 5 ] 

Fig-0 also displays as purple lines the sensi¬ 
tivities of the most relevant present (solid) and 
prospective 3000/fb (dashed) searches, namely 
those for t\ —$■ tyf in the CMSSM, NUHM1 and 
NUHM2 cases, and for t-\ —► bxf followed by 
xf —► Xi + soft particles in the pMSSMIO case. 
We see that the current search does not impact 
the CMSSM, NUHM1 or NUHM2. In the case of 
the pMSSMIO, the solid purple line is the current 
95% CL limit from the t\ —> bxf search in 
assuming a 100% branching ratio. This analysis 
is sensitive to stop topologies with b quarks in the 
final state where the decay products of the sub¬ 
sequent xf decay are undetected, and was used 
in [ 5 ] to constrain compressed stop spectra. How- 


45 
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■ CMSSM: best fit, la, 2 a 
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Figure 5. The (m~±,TO^.o) planes in the CMSSM (upper left), the NUHMl (upper right), the NUHM2 
(lower left) and the pMSSMIO (lower right). The red and blue solid lines are the Ax 2 = 2.30 and 
5.99 contours. The solid (dashed) orange lines are the current and projected 3000/fb 95% CL exclusion 
sensitivities for X 1 X 2 W/Z + I/t searches, the green dashed lines the projected 3000/fb 95% CL 
exclusion sensitivity for a X 1 X 2 W/h+J/r search (both from 21/, and the magenta dashed line 
is the projected 3000/fb 95% CL exclusion sensitivity for X 1 X 21 X 1 X 1 — > r, t —> 2,3 t's +1/t searches 
(from J/2 \). 


ever, it is not sensitive to decays involving on-shell 
W bosons or t\ —> ex/. We conclude that future 
searches have the potential to explore parts of the 
Ti coannihilation regions of the CMSSM, NUHMl 
and NUHM2, and of the Xi coannihilation region 


in the pMSSMIO casep~| but no DM channel can 
be fully explored by LHC searches. 


^We recall that the h and Z funnels in the pMSSMIO 
could in principle be explored by future searches for invis¬ 
ible h and Z decays. 
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■ CMSSM: best fit, 1 a, 2 a 
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Figure 6. The ( 171 ^, 171 ^ 0 ) planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 
(lower left) and the pMSSMIO (lower right). The red and blue solid lines are the Ay 2 = 2.30 and 5.99 
contours. The diagonal black dashed lines correspond to 0 and = mt + 0 . In each of 

the CMSSM, NUH M1 and NUHM2 panels, the solid purple line is the current 95 % C L limit from the 
ti —> tffl search in )3( and the dashed purple line is the 3000/fb projection from The solid and 

dashed purple lines for the pMSSMIO, obtained using m and 1)6 \ respectively, are the LHC Run 1 95% 
CL limit and the projected 3000/fb 95% CL exclusion sensitivity with 3000/fb for a t\ —>■ b\ 1 search, 
assuming a 100% branching ratio and rn^± — m^o = 5 GeV. 


4.6. The Heavy Higgs Bosons 

We now study the differences between the dom¬ 
inant DM mechanisms in the pMSSMIO and the 


other models in the (Ma, tan (3) planes shown in 
Fig. [7] In the case of the CMSSM, the regions al¬ 
lowed at the 95% CL and preferred at the 68% CL 
(blue and red contours, respectively) are generally 
at considerably larger Ma than the LHC bound 
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(shown as a solid purple line) 31 0 The fi 
coannihilation mechanism dominates in a region 
around M A ~ 2000 GeV for tan /3 < 40, the H/A 
funnel dominates for tan /3 ~ 50, and there is an 
intermediate hybrid region. On the other hand, 
Xi coannihilation dominates in an arc at larger 
M A . In the NUHM1, the hybrid and xf coanni¬ 
hilation regions are greatly expanded, and values 
of Ma closer to the LHC bound are allowed at 
the 95% CL. In the NUHM2, on the other hand, 
essentially all values of M A consistent with the 
LHC bound are allowed at the 95% CL, the Xi 
coannihilation mechanism dominates over most of 
the (Ma, tan/3) plane, leaving f\ coannihilation 
to dominate for M A > 2000 GeV and tan /3 < 30. 
Finally, we see that in the pMSSMIO xy coanni¬ 
hilation dominates the 68% CL region, that there 
is also a region at tan /3 < 20 where fi coanni¬ 
hilation may be important, and that there are 
intermediate uncoloured regions where neither of 
these mechanisms dominate. The LHC bound on 
Ma is again saturated for tan /3 < 50. 

We have also estimated (not shown) the 
prospective LHC 95% CL exclusion sensitivity in 
the H/A plane with 300/fb of data for the m™ ax 
scenario, scaling the current limit (using results 
from 147] and [48]), and comparing with the es¬ 
timated limits in 49], where good overall agree¬ 
ment was found. We estimate that M A < 2 TeV 
could be explored for tan/3 ~ 50, reducing to 
Ma /$ 1 TeV for tan/3 ~ 20. This would cover 
much of the H/A funnel and hybrid regions in 
the CMSSM, portions of these regions in the 
NUHM1, and parts of the xy coannihilation re¬ 
gions in the NUHM2 and the pMSSMIO, though 
not the regions of f\ coannihilation in these mod¬ 
els. Table |T] also summarizes the observability of 
the heavy Higgs bosons A/H in the different sce¬ 
narios considered. 

5. Direct Dark Matter Searches 

We now turn to the capabilities of direct DM 
search experiments to cast light on the various 


12 This line was calculated assuming the specific MSSM 
m™ scenario with MguSY = 1 TeV [32], and is used to 
give a rough impression of the location of the direct heavy 
Higgs-boson search bounds. 


DM mechanisms. Fig. kd displays the (m^cqcr® 1 ) 
planes for the CMSSM (upper left), the NUHM1 
(upper right), the NUHM2 (lower left) and the 
pMSSMIO (lower right), where er® 1 is the cross 
section for spin-independent scattering on a pro¬ 
ton. Our computation of u® 1 follows the proce¬ 
dure described in i3|, and we have once again 
adopted for the 7r-nucleon a term the value 
SttJV = 50±7 MeV. As previously, the Ay 2 = 2.30 
and 5.99 contours are shown as red and blue 


lines. The sensitivities of the XENON100 33 


and LUX 34 experiments are shown as green and 


black lines, respectively, and the prospective sen¬ 
sitivity of the LUX-Zepelin (LZ) experiment 35 
is shown as a purple line: the projected sensi¬ 


tivity of the XENON1T experiment 36 lies be¬ 


tween the current LUX bound and the future LZ 
sensitivity. Also shown, as a dashed orange line, 
is the neutrino ‘floor’, below which astrophysical 
neutrino backgrounds would dominate any DM 
signal 371 (yellow region). 

In the CMSSM case, we see that the current 
XENON100 and LUX data already put strong 
pressure on models where the focus-point or x/y 
coannihilation mechanism dominates. There are 
borderline regions that are formally excluded by 
the cr® 1 data considered in isolation, but become 
permitted at the 95% CL in a global fit includ¬ 
ing other observables, and also due to the uncer¬ 
tainties in the calculation of er® 1 that have been 
included in the evaluation of the global y 2 func¬ 
tion 13.. We also see that the xf coannihilation 
region and most of the H/A funnel region would 
be accessible to the planned LZ experiment. How¬ 
ever, much of the f\ coannihilation region lies 
below the LZ sensitivity, though it could be ac¬ 
cessible to a 20-tonne DM experiment such as 
Darwin 38 . However, this region can be cov¬ 


ered in the complementary direct searches at the 
LHC, as discussed previously. In the case of the 
NUHM1, the xy coannihilation and H/A funnel 
regions still lie largely within the LZ range, and 
the pure H/A funnel region is also within reach 
of LZ. However, the f\ coannihilation and hybrid 
regions extend below the reach of LZ and Dar¬ 
win, even below the neutrino ‘floor’. Again, this 
region could be partially covered by the comple¬ 
mentary LHC direct searches. Similar qualitative 
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Figure 7. The (Ma, tan/3) planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 
(lower left) and the pMSSMIO (lower right). The red and blue solid lines are the Ay 2 = 2.30 and 5.99 
contours, and the solid purple line is the current LHC 95% CL exclusion in the M™ ax scenario. 


conclusions apply to the NUHM2, with the addi¬ 
tional observation that the small t\ coannihilation 
regions lie below the LZ sensitivity and straddle 
the neutrino ‘floor’. 

Finally, we see that whereas the region of the 
pMSSMIO parameter space that is favoured at 
the 68% CL lies within reach of the LZ experi¬ 
ment, as is the case for much of the xf coanni¬ 
hilation region, there are models in the Xi and 


fi coannihilation regions, as well as in the h and 
Z funnels and uncoloured regions where none of 
these mechanisms dominate, in which cancella¬ 
tions in the spin-independent matrix element may 
drive ct® 1 below the neutrino ‘floor’. It should 
be kept in mind here (see the discussion in @)> 
that these very low values of ct® 1 are due to can¬ 
cellations 39 between different contributions to 
the matrix element for spin-independent scatter- 
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Figure 8. The (m^cpcr® 1 ) planes in the CMSSM (upper left), the NUHM1 (upper right), the NUHM2 
(lower left) and the pMSSMIO (lower right). The red and blue solid lines are the Ay 2 = 2.30 and 5.99 
contours, and the solid purple lines show the projected 95% exclusion sensitivity of the LUX-Zepelin 
(LZ) experiment \3f( . The green and black lines show the current sensitivities of the XENONIOO J33^j 
and LUX \3J)f experiments, respectively, and the dashed orange line shows the astrophysical neutrino 
‘floor’ \37j, below which astrophysical neutrino backgrounds dominate (yellow region). 


ing on protons. In [5] it was shown that similar 
cancellations hold when the cross section for spin- 
independent scattering on neutrons is considered, 
instead of the proton case shown in Fig. [8] 

Table [l] also summarizes the observability of 
DM particles in direct searches in the different 
scenarios considered. We see a degree of com¬ 


plementarity between the LHC and direct DM 
searches. 

We have focused in this article on the prospects 
for direct searches for DM scattering. A comple¬ 
mentary probe of the properties of supersymmet¬ 
ric DM is through indirect detection, searching 
for the traces of DM annihilation in the Galaxy. 
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A number of recent works have focussed on this. 
For example, 50 has demonstrated that Fermi- 
LAT satellite limits on 7 -ray emission in dwarf 
spheroidal galaxies 51 do not currently affect the 
parameter space of the pMSSM, although they 
may do so in the future □ Constraints from 
IceCube and the HESS telescope have been in¬ 
vestigated in 53 54 . IceCube limits 55 do not 
currently affect the pMSSMIO parameter space, 
while HESS bounds 56 are primarily on pure 
wino states, which must have masses greater than 
a TeV 57 58 in order to provide a thermal relic. 
Since the mass of the lightest neutralino is low in 
our models due to the incorporation of the (g— 2 )^ 
constraint, this does not affect our fits. Accord¬ 
ingly we do not include these data sets in this 
work, but we plan on implementing likelihoods 
for these searches in future analyses. 


6. Summary and Conclusions 

We have analyzed in this paper the mechanisms 
that play dominant roles in bringing the relic neu¬ 
tralino density into the range allowed by cosmol¬ 
ogy in the CMSSM, the NUHM1, the NUHM2 
and the pMSSMIO. We have delineated the re¬ 
gions of the parameter spaces of these models 
where dominant roles are played by fi, ti or \i 
coannihilation, or funnel regions where the neu¬ 
tralino annihilates rapidly via the heavy Higgs 
bosons H/A , and also regions where the neu¬ 
tralino has a significant Higgsino component. In 
the CMSSM, the NUHM1 and the NUHM2 we 
find that different mechanisms operate in differ¬ 
ent regions of the parameter spaces, with rela¬ 
tively small hybrid regions where two mechanisms 
contribute. In the pMSSMIO, xt coannihilation 
dominates in most of the parameter space, with 
some contributions from other processes in spe¬ 
cific ragions. 

Our assessments of the observability of super- 
symmetry within different models, depending on 
the dominant DM mechanisms, are summarized 
in Table [l] Within the CMSSM, the NUHM1 

13 Ref. |52| appeared while this work was being finalized, 
and may be used to constrain the Z- and fi-funnel regions 
of the pMSSMIO. We plan on incorporating this into a 
future analysis. 


and the NUHM2 ,^?t searches at the LHC can ex¬ 
plore significant portions of the fi coannihilation 
regions. These regions also offer the possibility 
that the fi may be relatively long-lived, and de¬ 
tectable at the LHC as a long-lived charged parti¬ 
cle. There are regions of the NUHM2 parameter 
space where t\ coannihilation dominates, which 
can also be explored by ]/Jt searches at the LHC. 
The X\ coannihilation and focus-point regions of 
these models can be explored by the LUX-Zepclin 
direct DM search experiment. Much of the H/A 
funnel regions in these models may be explored 
via LHC searches for the heavy Higgs bosons, and 
also via direct DM searches with the LUX-Zepelin 
and Darwin experiments. On the other hand, the 
fi coannihilation regions seem likely to lie beyond 
the reaches of the direct DM searches. 

Within our analysis of the pMSSMIO, xf 
coannihilation is the dominant DM mechanism 
in most of the parameter space, though this 
might change with a different set of independent 
pMSSM parameters. In addition, there are re¬ 
gions where fi coannihilation or direct-channel 
annihilation via a h and Z funnel may domi¬ 
nate. Parts of the pMSSMIO model space can 
be explored at the LHC via ]/Jt, H/A and other 
searches, and parts by direct DM searches. How¬ 
ever, we find no long-lived particle signature in 
the region of the pMSSMIO parameter space that 
is currently favoured statistically. Overall, large 
parts of the pMSSMIO parameter space could es¬ 
cape the LHC searches considered here, but a 
large fraction would be accessible to future DM 
experiments. 

Our analysis shows that the LHC and direct 
matter searches offer significant prospects for dis¬ 
covering SUSY if it is responsible for the cos¬ 
mological CDM, and in many cases the mode of 
discovery can reveal the nature of the dominant 
mechanism responsible for determining the CDM 
density. We look forward with interest to learning 
what the LHC and direct searches will be able to 
tell us about SUSY DM. 
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